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ABSTRACT 

We report on a short optical monitoring programme of the narrow-line Seyfert 1 Galaxy IRAS 
13224-3809. Previous X-ray observations of this object have shown persistent giant variabil- 
ity. The degree of variability at other wavelengths may then be used to constrain the conditions 
and emission processes within the nucleus. Optical variability is expected if the electron pop- 
ulation responsible for the soft X-ray emission is changing rapidly and Compton-upscattering 
infrared photons in the nucleus, or if the mechanism responsible for X-ray emission causes all 
the emission processes to vary together. We find that there is no significant optical variability 
with a firm upper limit of 2 per cent and conclude that the primary soft X-ray emission region 
produces little of the observed optical emission. The X-ray and optical emission regions must 
be physically distinct and any reprocessing of X-rays into the optical waveband occurs some 
distance from the nucleus. The lack of optical variability indicates that the energy density of 
infrared radiation in the nucleus is at most equal to that of the ultraviolet radiation since lit- 
tle is upscattered into the optical waveband. The extremely large X-ray variability of IRAS 
13224-3809 may be explained by relativistic boosting of more modest variations. Although 
such boosting enhances X-ray variability over optical variability, this only partially explains 
the lack of optical variability. 
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1 INTRODUCTION 

X-ray observations by Boiler et al (1997) have shown persistent 
giant variability on short timescales, with an amplitude of vari- 
ability far in excess of that seen in a typical broad line Seyfert 
1 galaxy. Such extreme variability is determined by the emission 
and variability mechanisms in the nucleus. The degree of variabil- 
ity at other wavelengths may be used to constrain the conditions 
and mechanisms. Optical variability will occur if, for example, the 
X-ray emitting electron population is rapidly changing and Comp- 
ton scattering infrared radiation in the nucleus, or if the mecha- 
nism responsible for the X-ray variations causes all the emission 
processes to vary together. It has also been suggested (Boiler et al 
1997) that such large amplitude X-ray variability may be due to 
relativistic boosting, the degree of which is spectrally dependent. 
We have therefore observed IRAS 13224-3809 on three consecu- 
tive nights in order to determine the level of optical variability in 
the source. 

Previous studies of the narrow line Seyfert 1 galaxy IRAS 
13224-3809 have shown it to be variable both in the UV and X- 
ray wavebands. Long timescale Lya line variability has been ob- 
served, with the line profile and flux changing between three IUE 
observations over four months (Mas-Hesse et al. 1994). UV contin- 



uum variability has been observed in 1 1 IUE observations spanning 
three years, with the flux varying by 24 per cent (Rodriguez-Pascual 
1997). 



2 OPTICAL OBSERVATIONS 

IRAS 13224-3809 (q 20 oo = 13 h 25 m 19 s , 5 2 ooo = -38 h 24 m 53 s ) 
was observed in B, V, R and / on three consecutive nights from 
1997 March 13 to 16 using the blue sensitive TEK4 CCD at the 
Cassegrain focus of the lm lacobus Kapteyn Telescope (JKT). The 
5.6 x 5.6 arcmin CCD field contained the galaxy and a number 
of foreground stars (see Fig. 1) allowing us to perform relative 
photometry. The images were bias subtracted and flat fielded us- 
ing sky flats. Flux calibration was performed using the standard 
stars PG0942-029B, C and D of Landolt (1992). IRAS 13224-3809 
has a foreground star 7.6 arcsec from the optical nucleus (10.7 arc- 
sec from the X-ray centroid position (Boiler et al. 1997)), and the 
aperture used to collect the light from the nucleus was chosen to 
have a diameter of 6.6 arcsec. This was found to collect most of the 
light from the nucleus whilst minimising that from the surrounding 
galaxy and the foreground star. Repeating our analysis with differ- 
ent sized apertures did not significantly affect our results. The same 
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Figure 1. V-band image of IRAS 13224-3809 (a 2 ooo = 13 h 25 m 19 s , 
52000 = — 38 h 24 m 53 s ) and a comparison star A. 
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Figure 2. Radial profile of IRAS 13224-3809. The effect of the point source 
close to the nucleus can be seen at a distance of about 23 pixels from the 
centre. The scatter in pixel value is due to the contribution from different 
azimuths. 



size aperture was used to obtain the flux from a comparison star in 
the field. Fig. 2 shows the radial profile of pixel values from the 
centre of IRAS 13224-3809. The scatter in pixel value is due to the 
different contributions to the radial profile from different azimuths. 
The contribution from the nucleus, underlying galaxy, the nearby 
star and background may be seen. 

We follow the method of Done et al (1990) in which the vari- 
ations seen in a comparison star due to small-scale rapid changes 
in seeing are used as a scalable template to model the variability 
one would expect to see in the nucleus as a results of those see- 
ing changes. During each night the zenith distance of IRAS 13224- 
3809 varied between 67-72° and over such a small range we expect 
the extinction to remain almost constant. Since the images were 
taken through a number of different filters we grouped the images 
by night and by filter and then normalised them so that they have 
the same mean. This allows us to compare images taken with dif- 
ferent filters. Fig. 3 shows the stellar and galactic residuals about 
the mean for a typical night. We would expect the slope of the best 
fitting line to be 1 if the seeing had an equal effect on both the 
galaxy and the comparison star. The actual fit has a slope of 1.10 
indicating that seeing has a larger effect on the flux measured for 
the galaxy than the star, assuming both are not variable. The stellar 
residuals are then used as a scalable template for the variations of 
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Figure 3. The residuals about the mean for the nucleus and star for the 
second night's data. The best fitting straight line has a slope of 1.10. This 
indicates that (were both sources to be constant) the seeing has a larger 
effect on the galaxy than the star. 
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Figure 4. Flux of the nucleus predicted using the scaled stellar residuals 
(solid line) and the actual observed flux of the galaxy. 

the galaxy about the mean flux of the galaxy. The normalisation of 
this template is chosen to minimise the difference between the pre- 
dicted and observed galactic fluxes. Fig. 4 shows the template light 
curve, created using the scaled stellar residuals, and the actual light 
curve. There is very good agreement between the two. 

Fig. 5 shows the light curve for the three nights using data 
from all filters, with la Poisson error bars. The light curve is con- 
sistent with less than 1 per cent variability during each night. One 
may expect the amplitude of variability to be wavelength depen- 
dent, but we are unable to detect significant variability in the data 
from any one filter. The light curve for each night is normalised to 
the mean for that night although there is less than 1 per cent vari- 
ability between nights. The large breaks in the light curve are due 
to high winds at the telecsope preventing observing. 



3 DISCUSSION 
3.1 Optical 

The optical light curve of IRAS 13224-3809 is consistent with the 
short time-scale variability of the source being less than 1 per cent. 
We are unable to comment on the wavelength dependence of any 
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Figure 5. The light curve for the three nights. It is consistent with the source 
having varied by less than 1 per cent during the course of our observations. 



optical variability. Our aperture, however, also collects flux from 
the underlying galaxy which dilutes the signal from the nucleus. 
As may be seen from Fig. 1 the three main components of the flux 
received within a radius of 10 pixels, namely the contribution of 
the background, the underlying galaxy and the nucleus may be sep- 
arated. If the nucleus is assumed to be a Gaussian point source and 
we (over-) estimate the contribution from the galaxy then the vari- 
ability of the nucleus may be constrained. It is found that a 1 per 
cent variation of the entire source may correspond to a ~ 2 per cent 
variation in the nucleus. 

The continuum optical spectrum of IRAS 13224-3809 (Boiler 
et al 1993) has a photon index of approximately r o ~ across the 
optical waveband. The V band magnitude of the nucleus is 15.2, 
corresponding to a vF v flux in the V band of 1.5 x 10~ n erg 
cm~ 2 s" 1 at5500A. 
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Figure 6. The ROSAT 30-day HRI light curve of IRAS 13224-3809 (Boiler 
et al 1997) in linear scale (upper panel) and logarithmic scale (lower panel). 
Note that the source varies cosiderably even during apparently quiet be- 
haviour in the upper panel, e.g. days 8-10. 



3.2 X-ray 

Simultaneous X-ray data are not available but if the 30 day ROSAT 
light curve of Boiler et al (1997) is assumed to be typical of IRAS 
13224-3809 it is extremely unlikely that the X-ray flux did not at 
least double during the period of our observations. This light curve 
is shown in Fig. 6 and there is no set of three consecutive days dur- 
ing which the X-ray flux did not change considerably. The only pe- 
riod of relative 'quiescence' is between 8-10 days, and even during 
this the X-ray flux at least doubles. A more recent HRI monitoring 
campaign has confirmed the continued existence of such extreme 
variability (private communication). 

The X-ray photon index of IRAS 13224-3809 was observed 
to be very high, T x ~ 4, and the average HRI count rate was 4.7 x 
10~ 3 count s _1 (Boiler et al 1997) corresponding to a vF v flux of 
2.9 x 10~ 13 erg cm" 2 s" 1 at 1.25 keV. 

3.3 The nuclear X-ray and optical emitting regions 

Optical emission can originate from the soft X-ray emission re- 
gion, or very close to it, by several means. It may be emitted by the 
same electrons that produce the X-rays, by cyclotron or synchrotron 
emission or Comptonization (perhaps of photons produced by the 
cyclo-synchrotron process; see e.g. Di Matteo, Celotti & Fabian 
1997). Comptonization is however unlikely to make a large contri- 
bution to the observed flux since the steep X-ray spectrum implies 



a low Compton {/-parameter which gives a similarly steep optical 
spectrum. Optical emission may also originate from reprocessing of 
the X-ray flux, and therefore would have a thermal spectrum (prob- 
ably from a range of temperatures). Finally it may be triggered by 
whatever mechanism causes the X-ray variability but not share the 
same emitting electrons as the X-rays. 

The X-ray flux is expected to have at least varied by a factor 
of two during our observations. Let us assume that the X-ray flux 
doubled at some point and consider the implications of the lack of 
optical variability. 

A lack of variability in the optical band can be explained in 
several ways. It may simply be due to most, more than 98 per cent, 
of the optical emission emerges from a region unconnected with 
that producing the soft X-ray flux. Note however that the average 
soft X-ray flux (below 1 keV) is comparable to the optical flux and 
the peak soft X-ray flux may be significantly greater. It could also 
be due to the presence of many very dense clouds in the nucleus 
which free-free absorb the intrinsic optical flux (Celotti, Rees & 
Fabian 1994). Alternatively, the spectrum of the X-ray source may 
be such that relativistic effects preferentially enhance the variability 
of the X-ray emission over the optical emission. 

The exceptionally large and rapid X-ray variability of IRAS 
13224-3809 may be explained if the X-ray emission is relativisti- 
cally boosted (Boiler et al 1997). For a source moving at a frac- 
tion 13 of the speed of light inclined at an angle i to the observer 
the Doppler parameter is given by S = [7(1 — /3sini)] _1 , where 
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Figure 7. X-ray (T x = 4, solid) and optical (T = 0, dashed) boost factors 
as a function of radius. Curves are shown for two disc inclinations, 30° 
(lower curves) and 90° (upper curves). Velocities in the disc were calculated 
using a Newtonian pseudo-potential. 



7 = (1 — f3 2 )~ If the source is a power law of photon index T 
this gives rise to a boost in the amplitude of variability by a fac- 
tor of S i+F . Fig. 7 shows the X-ray and optical boost factors that 
may be produced as a function of radius for different continuum 
photon indices and accretion disc inclinations. It is possible for the 
X-ray boost factor to be many times greater than the optical boost 
factor. (The ratio of the X-ray to optical boost factor is ~ 5 rx ~ r °). 
The fraction of the optical emission that may be produced by, or 
whose variability is tied to, the soft X-ray emitting regions may 
then be constrained. The maximum fraction 0.025 rx ~ r ° of the op- 
tical emission that is due to the X-ray source is shown in Fig. 8 
assuming values of T x —4 and r o = 0. For reasonable values 
of the inclination and emission radius it is unlikely that the X-ray 
source is responsible for more than about 10 per cent of the optical 
flux. 

The observed photon index of the optical continuum may dif- 
fer from that of the optical emission associated with the X-ray 
source. The lowest possible T is — 1 from the Rayleigh- Jeans part 
of the blackbody spectrum, and it is possible that T x exceeds 4. It 
remains unlikely, however, that more than about 20 per cent of the 
optical emission is closely associated with the soft X-ray emission. 
Differential relativistic effects between the X-ray and optical bands 
can therefore be up to a factor of 10 but not completely explain the 
lack of optical variability. 

The hot electron population responsible for the X-ray emis- 
sion can inverse-Compton-scatter lower energy infrared photons 
into the optical waveband and similarly ultraviolet photons into the 
X-ray waveband. The observed lack of variability can be used to 
place a limit on the ratio of the energy density of infrared to ultra- 
violet radiation in the nucleus. The observed ratio of 1.25 keV X- 
ray to optical luminosity, L x /L ~ 10~ 2 , and the lack of optical 
variability suggests that the variable X-rays have Z/ X ( va r) /£ (var) 
less than or comparable with 1, so, in the nucleus, the energy den- 
sity of infrared radiation £ir is at most that of the ultraviolet ra- 
diation £tjv- If we assume the source to have an X-ray photon 
index T x = 4, the 1.25 keV flux implies the 0.1 keV vF v flux 
will be 4.5 x 10 -11 erg cm~ 2 s _1 allowing the stronger constraint 
em < 0.01 x £uv 
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Figure 8. For a given radius and inclination the maximum fraction of the 
optical emission that may be produced by the X-ray source is calculated. 
The solid curves are for inclinations of 30° , 60° and 90° (increasing up- 
wards). This assumes that the X-ray flux is seen to double and the optical 
flux increases by at most 2 per cent. The dotted curve shows the limit on the 
optical fraction if the variability is not due to relativistic effects. 



3.4 Location of any reprocessing material 

The lack of optical variability implies that very little of the primary 
soft X-ray flux is reprocessed into the optical waveband within ap- 
proximately a light day of the source. This may be due to the in- 
ner disc where the X-rays originate being hot and small, radiating 
mostly in the EUV, and to a lack of optically thick material oriented 
to reprocess these X-rays at larger radii. For a 3 x 1O 7 M0 black 
hole two light days corresponds to roughly 1000 Schwarzschild 
radii. 



4 CONCLUSION 

The Narrow Line Seyfert 1 galaxy IRAS 13224-3809 shows persis- 
tent giant amplitude X-ray variability yet almost no optical variabil- 
ity exceeding 2 per cent. Such an extreme difference in variability 
between the two wavebands suggests the X-ray emitting regions do 
not produce any optical variability. This is expected if the electron 
populations responsible for X-ray and optical emission are phys- 
ically distinct and the spectrum of the X-ray source intrinsically 
produces almost no optical emission. The energy density of infrared 
radiation in the nucleus is at most equal that of the ultraviolet ra- 
diation since the X-ray emitting electrons are not responsible for 
significant inverse Compton scattering into the optical waveband. 
If relativistic boosting occurs then the conclusions are weakened 
somewhat with at most about 20 per cent of the optical emission as- 
sociated with the X-ray source. Any significant reprocessing of pri- 
mary soft X-rays into the optical waveband occurs on scales greater 
than 1000-Rs The lack of optical variability is consistent with ob- 
servations of less extremely X-ray variable Seyfert 1 galaxies such 
as NGC 4051 (Done et al 1990). 
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